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We investigate the thermalization and the chemical equilibration of a parton plasma created from 
Au+Au collision at LHC and RHIC energies starting from the early moment when the particle 
momentum distributions in the central region become for the first time isotropic due to longitudinal 
no . cooling. Using the relaxation time approximation for the collision terms in the Boltzmann equations 

on ■ 

for gluons and for quarks and the real collision terms constructed from the simplest QCD interactions, 
we show that the collision times have the right behaviour for equilibration. The magnitude of the 
quark (antiquark) collision time remains bigger than the gluon collision time throughout the lifetime 
Oh' °f the plasma so that gluons are equilibrating faster than quarks both chemically and kinetically. 

That is we have a two-stage equilibration scenario as has been pointed out already by Shuryak 
sometimes ago. Full kinetic equilibration is however slow and chemical equilibration cannot be 
completed before the onset of the deconfinement phase transition assumed to be at T c = 200 
MeV. By comparing the collision entropy density rates of the different processes, we show explicitly 
that inelastic processes, and not elastic processes as is commonly assumed, are dominant in the 
equilibration of the plasma and that gluon branching leads the other processes in entropy generation. 
We also show that, within perturbative QCD, processes with higher power in a s need not be less 
important for the purpose of equilibration than those with lower power. The state of equilibration of 
the system has also a role to play. We compare our results with those of the parton cascade model. 
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A goal of the future heavy ion collision experiments at the relativistic heavy ion collider (RHIC) at Brookhaven 
and at the large hadron collider (LHC) at CERN is to find the quark-gluon plasma. The primary aim is of course to 
show that quarks and gluons can indeed be freed from their hadronic "prison" and exist as individual entities in a hot 
plasma. Once this is realized, one can then turn to the diverse physics of such a new state of matter. One of these is 
the relation of the various thermodynamic variables to each other or in other words, the equation of state [JO. In order 
to probe this in experiments, an equilibrated quark-gluon plasma is required. In this work, we look at how far can one 
expect to have such a plasma in equilibrium. Because of the importance of this question, various different approaches 
have already been taken to address this issue. In particular, Shuryak ||] argued that equilibration of the plasma 
proceeds via two stages in the "hot gluon scenario" . First the equilibration of the gluons and then that of the quarks 
follows with a certain time delay. Thermal equilibration is quite short for gluon < 1 fm with high initial temperature 
of 440 MeV at LHC and 340 MeV at RHIC. However, these estimates are based on thermal reaction rates for large 
and small angle scatterings and on the assumption that one scattering is sufficient to achieve isotropy of momentum 
distribution. As has been shown in || using a family of different power behaviours for the time-dependence of the 
collision time, the assumption of one scattering is sufficient is a serious underestimate. With a larger number of 
scatterings, using the same arguments as in the initial temperature will be lowered and the thermalization time 
will be increased. Also, we argue that estimates based on using the scattering rate alone is incorrect, since in a 
medium, one must consider the difference of the scattering going forward and backward both weighed with suitable 
factors of particle distribution functions. Hence, the process with the largest cross-section is not necessarily the more 
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important. However, we will show the two-stage equilibration scenario or in other words, gluons equilibrate much 
faster than quarks and antiquarks. 

The other approach is the semi-classical parton cascade model (PCM) 0-0], which is based on solving a set 
of relativistic transport equations in full six-dimensional phase space using perturbative QCD calculation for the 
interactions, predicts an equilibration time of 2.4 fm/c for Au+Au collision at 200 GeV/nucleon. This approach, 
which uses a spatial and momentum distribution obtained from the measured nuclear structure functions for the 
partons as initial state, is very complicated. Due to the finite size of the colliding nuclei, it is hard to clearly identify 
thermalization in terms of the expected time-dependent behaviours of the various collective variables . But by fitting 
the total particle rapidity and transverse momentum distributions of the defined central volume, roughly identical 
temperatures are obtained || and hence the claim of thermalization. However, in terms of the same distributions of 
the individual parton components, this becomes less obvious to be the case J(|. As was stated in Q, the momentum 
distributions are not perfect exponentials and therefore there is no complete thermalization in any case. 

We will look at this problem of equilibration using a much simpler approach which is based on the Boltzmann 
equation and the relaxation time approximation for the collision terms. Initially used by Baym to study thermal 
equilibration and has subsequently been used in the study of various related problems J8|-[ll| . The conclusion of these 
works is, in general, if the collision time 6 which enters in the relaxation approximation, grows less fast than the 
expansion time r, then thermal equilibration can be achieved eventually. In the case of the quark-gluon plasma, it is 
not sufficient to know that equilibration will be achieved eventually because the plasma has not an infinite lifetime in 
which to equilibrate. We would like to know how far can it equilibrate before the phase transition. To answer such a 
question, we will use both the relaxation time approximation and the interactions obtained from perturbative QCD 
for the collision terms to determine 0. This approach has been used previously to study both thermal and chemical 
equilibration in a gluon plasma [ ]TT[ | where it was found that with the initial conditions obtained from HIJING results, 
the gluon plasma had not quite enough time to completely equilibrate. In the present case of a quark and gluon 
parton plasma, quarks and gluons are treated as different particle species rather than as generic partons and so they 
have different time-dependent collision times. As a result, they approach equilibrium at different rates and towards 
different target temperatures. The latters will converge only at large times. It follows that the system can only 
equilibrate as one single system at large times. This lends support to the two-stage equilibration scenario B. 

In an expanding system, particles are not in equilibrium early on because interactions are not fast enough to 
maintain this so they are most likely to start off free streaming in the beam direction Jl0| , p^[ . Thermalization will be 
seen as the gradual reduction of this free streaming effect as interactions gain pace and momentum transfer processes 
are put into action to bring the particle momenta into an isotropic distribution. The present approach takes into 
account of these effects. 

As in the previous work pjj ], isotropic momentaneously thermalized initial conditions are used at both RHIC 
and LHC energies. These are obtained from HIJING results after allowing the partons to free stream until the 
momentum distribution becomes isotropic for the first time |T^-|T^]. From then on, interactions are turned on but 
the distribution becomes anisotropic again due to the tendency of the particles to continue to free stream. It is the 
role of interactions to reduce this and to progressively bring the distributions into the equilibrium forms. We have 
shown that, surprisingly, kinetic equilibration in a pure gluon plasma is driven mainly by gluon multiplication and 
not gluon-gluon elastic scattering. In this paper, we include quarks and antiquarks and consider the equilibration of 
a proper QCD plasma. We explicitly break down the equilibration process into each of its contributing elements and 
show which interactions are more important and hence uncover the dominant processes for equilibration. In fact, our 
result is inelastic interactions are most important for this purpose both for quarks and for gluons. 

Our paper is organized as follows. In Sect. O , we describe the Boltzmann equations with the relaxation time 
approximation for two particle species. In Sect. |III| the time-dependent behaviour of the collision times, 0's, necessary 
for equilibration will be analysed and extracted. The particle interactions entering into the collision terms and details 
of their calculations will be explained in Sect. IV. Initial conditions used will be given in Sect. |y| and lastly the 
results of the evolution of the plasma will be shown and discussed in Sect, [vl We finish with a brief discussion of 
the differences with the results of PCM. 



II. RELAXATION TIME APPROXIMATION FOR TWO PARTICLE SPECIES 

In the absence of relativistic quantum transport theory derived from first principle of QCD p"^-|22|], we base our 
approach on Boltzmann equation with both the relaxation time approximation for the collision terms and the real 
collision terms obtained from perturbative QCD. Treating quarks and gluons on different footings, we write down the 
Boltzmann equations 
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(1) 



where fi is the one-particle distribution and stands for the collision terms and includes all the relevant interactions 
for particle species i and i = g,q,q. Concentrating in the central region of the collision where we assumed to be 
spatially homogeneous, baryon free and boost invariant in the z-direction (beam direction) so that f q — f q and 
fi = fi(p±,p' z ,r) where p' z = ^{p z — up) with 7 = 1/yl — u 2 and u = z/t is the boosted particle z-momcntum 
component and r = y/t 2 — z 2 is the proper time. Following Baym 0, the Boltzmann equation can be rewritten as 



dfi 
dr 



= C l (p±,P z ,r) (2) 

p z T 



in the central region. Using the relaxation time approximation 

Ci{p±,p z ,T) = — — (3 

9i{T) 

where f eq i is the equilibrium distribution and Oi is the collision time for species i, this allows us to write down a 
solution to Eq. (|J). 

Mp,t) = f 0i (p±,p*T/T )e- Xi + £ dx\e<-^! eql ^p\ + (p z rlT') 2 ,T eql {r')) , (4) 

where 



f 0i (p±,p z T/r ) = (eM^Pl + (Pzt/t ) 2 /To)/1 t T l) , (5) 

is the solution to Eq. (Q) when C — which is also the distribution function at the initial isotropic time To , with initial 
fugacities Iq and temperature T . It is of such a form because of the assumption of momentaneously thermalized 
initial condition. The functions a;j(r)'s, given by 



i(r) = f dr'/diir') 

J Tn 



(0) 



play the same role as O^s in the sense that their time-dependent behaviours control thcrmalization. T eq that appears 
in f eq i, is the time-dependent momentaneous target equilibrium temperature for the i particle species. The two terms 
of equation Eq. (|J) can be thought of, up to exponential factor, as the free streaming (first term) and equilibrium 
term (second term). Whether species i equilibrates or not depends on which of the two terms dominates. 

In the present case of two species, the energy conservation equations are, in terms of the equilibrium ideal gas 
energy densities e eq g — aiT^ q g , e eq q = nfb2T^ q q , a-i = 87r 2 /15, 62 = 77r 2 /40 and rif is the number of quark flavours, 

d^-l ^ £j ~\~ PL i ^eq i 

dr T 6i 

and 

detot _^ (tot + Pl tot 
dr t 

where e tot = J2i e * an d PL tot = HiPh i, or in other words 

e; - e e 



(8) 



Et^ = °- ( 9 ) 



The above equation only expresses the fact that energy loss of one species must be the gain of the other. The transport 
equations of the different particle species are therefore coupled as they should be. The longitudinal and transverse 
pressures are defined as before 

PL,T i{r)=Vi [ — fi(p±,Pz,r) , (10) 

J (27T) d p 
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with v g = 2 x 8 = 16 and i/ ? = 2 x 3 x nj = 6nj, the multiplicities of gluons and quarks respectively. 

Here the equilibrium target temperatures T eq g and T eq q cannot be the same in general since, as we will see in 



Sect. VI, 9 g ^ 8 q = 9q. Therefore gluons and quarks will approach equilibrium at different rates. Note that energy 



conservation here does not mean 

£g ~t~ € e q g -f- 2e e q q (1 1 ) 

since 9 g < 9 q always, at least at small times, so gluon energy density e g will approach e eq g faster than e q approaches 
e eq q so the two equilibrium energy densities should not be considered to be those which can coexist at the same 
moment. This can only be true at large t when T eq g ~ T eq q and 9 g ~ 9 q . If Eq. ( pi] ) were true, the condition 
for energy conservation Eq. (||) could not hold when 9 g ^ 9 q . Since our QCD plasma is a dynamical system under 
one-dimensional expansion as well as particle production, the target temperatures T eq g and T eq q must be changing 
continuously and must approach each other at large times before the gluon and quark (antiquark) subsystems can 
merge into one system and exist at one single temperature. Likewise, we believe 9 g and 9 q should also converge to a 
single value at large times, unfortunately, this will take too long to happen in the evolution of our plasma although 
we can be sure that both 9 g and 9 q increase less fast than the expansion time t near the end of the evolution, a 



condition which, as has already been stated in the introduction and we will see again in Sect, til, is necessary for 
thcrmalization. 



Before considering the evolution of the QCD plasma under real interactions, we can deduce analytically, using 
Eq. (2) and Eq. (^), the conditions on the 0j's under which the plasma will come to kinetic equilibrium. Multiplying 
Eq. (4) by particle energy and integrating over momentum, we have the equations for the e^'s. Further manipulating 
these gives, 



dx 



o 



< e<{r'h{r' /r)(e eqi {r') - e,(r')) - ^ (t'^t'/t)^')) } = , (12) 



where 



sin 1 \J\ — r2 



h(r)= / dyy/l-y*(l-r*) = -[r+- 7 =L= (13) 



o 2 V VT^ 1 

and x[ — Xi(r'). Supposing as r — > oo, x g — > oo and x q — > oo then the integrand in Eq. ( |i~2] ) will be weighed by 
the r' — > oo or large x\ limit. It follows that the term within braces in Eq. ( |l2|) must be zero at large r' so using 
h'{r)\ r= i — 1/3, we have 

d^i 4 e eq i . . 

This means each species will undergo near hydrodynamic expansion at large r modified by energy lost to or energy 
gained from the other species. The latter should be small at such times. Summing Eq. ( p^ ) over species, we obtain 
the energy conservation equation for a system undergoing hydrodynamic expansion 

^i + |^=0, (15) 

GtT 3 T 

with p L tot = e tot /3. 

If one 9i is such that the corresponding Xi — > xi x < oo as r — > oo then hydrodynamic expansion does not apply to 
that species since we have 

d(eir) _ (ej - e eq i) 



Phi , (16) 



dT 6i 

where now pl i ^ so kinetic equilibrium is not established. The r.h.s. of Eq. ( [l6| ) is negative if these particles are 
losing energy or gaining energy at a rate less than i/r at large r. Therefore e^r must decrease towards a non-zero 
asymptotic value (e^r),^, since Xi oo < oo e;r > always, which results in a free streaming final state for these 
particles 
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€i(r -> oo) ~ (e i r) 00 /r . (17) 

A similar free streaming final state will be reached if the rate of gaining energy is larger than pl %/t at large r. In 
this case, although e^r is increasing, the 6j of the other particle species with Xj — * oo as r — ► oo will be close to e eq j 
and so the energy transfer will be very small. One can deduce that as r — > oo 

l»££i^il^ >^ =► ^^0, (18) 

hence e^r — * (eir)oo- That is ejT now increases towards some asymptotic value instead of decreasing towards one as 
in the previous case. But it ends up with a free streaming final state nevertheless. We do not consider the case where 
the relative rate (e eq j — e^r/Oip^ $ oscillates about one at large r except to say that on the average d(ejT)/dr ~ 
and so an average free streaming final state is likely. 

The last possibility where x, — > oo < oo as r — > oo for both particle species, Eq. ( |l6| ) applies to both. Barring 
the case of the oscillating relative rate, one particle species must lose energy and so by the above argument, a free 
streaming final state results. For the remaining particle species, it does not matter whether d(eiT)/dr is or is not 
positive at large r, these particles will also be in a free streaming final state. If the rate is negative, then the same 
argument that leads to Eq. ( |l7| ) applies. If it is positive, since the species that is losing energy is approaching free 
streaming so the energy transfer must go to zero. Then we are back to Eq. (|l8|). 

The conclusions are therefore, depending on the time-dependent behaviours of 9 g and 9 q , 

1. x g — > oo and x q — ► oo as r — * oo are required for the whole system to completely thermalize. 

2. x g — > oo and x q — > x q oo < oo or x q — > oo and x g — > x g oo < oo as r — > oo imply that only the species with 
Xj — > oo will thermalize, the other species will not equilibrate but free streams at the end. The system will end 
up somewhere between free streaming and hydrodynamic expansion. 

3. Both x g — > x g oo < oo and x q — > x q oo < oo as r — > oo then the whole system will end up in a free streaming 
final state. 

One can understand these Xi behaviours in terms of 9^s by assuming simple power r-dependence for the latters. 
One finds that 9^s must all grow slower than r for the whole system to achieve thermalization. If cither one or more 
grow faster then a mixed or a complete free streaming final state results. 

IV. PARTICLE INTERACTIONS — COLLISION TERMS 

To investigate the evolution of a proper QCD plasma, we consider the following simplest interactions at the tree 
level 

99 < — * 999 , 99 < — » 99 , (19) 

99 < — > qq , gq < — > gq > gq < — » gq > (20) 

qq < ► qq , qq < ► qq , qq < > qq . (21) 

As in [^3|-[T^] , we include only the leading inelastic processes i.e. the first interaction of Eq. ( |l9| ) and Eq. (pfj|)p1. We 



will return to this point later on in Sect. VI 



In the solutions Eq. (|4|) to the Boltzmann equations Eq. (|2j), there are two time-dependent unknown parameters 
9i and T eq , for each species which very much control the particle distributions. To determine them, we need two 
equations each for gluons and for quarks. In order to show the relative importance of the various interactions Eq. (|lS|), 
( po| ) and (^l[) in equilibration, we find these time-dependent parameters by constructing equations from the rates of 
energy density transfer between quarks (antiquarks) and gluons and the collision entropy density rates. 



x The first one of Eq. ( |2l| ) could also be inelastic but here we give the same chemical potential to all the fermions so we do not 
consider quark-antiquark annihilations into different flavours as inelastic for our purpose. 
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From Eq. (§), (§) and @, the energy density transfer rates are 



dr t 6i J (2tt) 



pCi( P± , Pz ,T)=£ z , (22) 



where £i is the energy gain or loss of species i per unit time per unit volume. As stated in Sect. O, Si's must obey 
J2i Si — for energy conservation. 

The other equations, the collision entropy rates can be deduced from the explicit expression of the entropy density 
in terms of particle distribution function p3fl 

*(r) = -i/i J ^{/ l (p,T)ln/ 4 (p,r)T(l±/ i (p,r))ln(l±/ J (p,r))} , (23) 

where the different signs are for bosons and fermions respectively. They are, using again Eq. (0), (||) and (0), 

(S 1 ) co « = ~ Vl !wf (^L ^ ( IT^) (24) 

wfNf 11 ^^)- (26) 



By using the explicit expression for the collision terms Cj's constructed from the interactions Eq. (|19|), ( |20| ) and 
( pl| ) within perturbative QCD, Eq. (p2|), ( p5| ) and ( prjj ) allow us to solve for 0,'s and T eg j's. 

The gluon multiplication contribution to C g is constructed from the infrared regularized Bertsch and Gunion 
formula [^4| for the amplitude with partial incorporation of Landau-Pomeranchuk-Migdal suppression (LPM) for 
gluon emission and absorption fl3| , ^5| -p7| as in the previous work |TT1 . The explicit form of the gluon multiplication 
collision term and a discussion of the problem regarding how to incorporate the LPM effect correctly can be found 
there also. The remaining binary interaction contributions to Ci for particle 1 is, as usual, given by 

x[/i/ 2 (l ± / a )(l ± h) - hh{\ ± A)(l ± h)\ (27) 

where the Vi runs over all the binary processes in Eq. (JTs|) , ( ]20| ) and ( pi] ) which involve species i, {Ai^l 2 is the sum 
over final states and averaged over initial state squared matrix element, S-p i is a symmetry factor for any identical 
particles in the final states for the process Vi and v 2 is the multiplicity of particle 2. 

We take \M. Vi \ 2 's from [^8| and infrared regularized them using either the Debye mass m 2 D for gluons or the quark 
medium mass m 2 . for quarks to cut off any infrared divergence. These masses are now time-dependent quantities in a 
non-equilibrium environment. With non-isotropic momentum distribution, both the Debye mass p9| , p^| and the gluon 
medium mass, m 2 , are directional dependent. This is, however, not the case for the quark medium mass, m 2 , which 
remains directional independent as in equilibrium. The directional dependence arises out of the cancellations between 
identical type of distribution functions similar to those one finds in the derivation of hard thermal loops 1 3 1 32] |. To 
keep things simple, we removed the directional dependence from m 2 D and use, for SU(N=3), to leading order in a s , 



f d^ d / 
m%(r)=-8na s J -^ — (n f g + n f f q ) 



(28) 



For the quark medium mass, to the same order, we use 

'N 2 - 1\ f d 3 p 1 



*>**(-iF) < 2<J) 



(2»)» |p 

which is just the equilibrium expression but with non-equilibrium distribution functions. 

With these masses, we regularize the squared matrix elements by hand and inserting the masses as follows. 



G 



, . l9 9 q 2 ( ut us st , 



2 V (s + ra|,) 2 (t-m 2 D ) 2 (u-m 2 D ) 2 

\ M \*- 9 l( —t | u ^1 3 " 2+ * 2 (31) 

i^™i 2 = -i 2 - 52 ( x - (t^b* " * ( vf^) + i^i) ) ) (32) 

2 ... ,2 2 5 2 /2( S 2 + t 2 ) , x 2( M 2 + S 2 ) 



|M g ^ M | - \M qq ^ M \ - — ( v(u _ to 2 ))2 + Su (t _ m 2 D y 



As 2 

~~ o 71 2"77 2~7 (^3) 

2 3 2 / s 2(s 2 + t 2 ) , . x 2(t 2 + « 2 ) 



l^™! 2 " ^ ^ 13 ^ 24 ^# + ' 5l ^ 34 ( S + 4m 2 ) 2 



4 i 2 



-^3*34 3 {u _ ml)[s + Aml) ) ( 34 ) 

where the Sij signifies that the i and j quark or antiquark must be of the same flavour. This regularization amounts 
to screening spacelike and timelike infrared gluons by m 2 D and 4m 2 , respectively and infrared quarks by m 2 . We stress 
that this regularization is done in a very simple manner and with the right order of magnitude for the cutoffs. Its aim 
is to get some estimates to the collision rates without involving too much with the exact and necessarily complicated 
momentum dependent form of the true infrared screening self-energies in an out-of-equilibrium plasma when their 
infrared screening effects should be in action. They should be the extension of the 2-point gluon and quark hard 
thermal loops |5l]-[35| to a non-thermalized environment. 

We should mention here that the choice of the pair of equations for solving the two time-dependent unknowns 0i 
and T eq , for each particle species is not unique. One can equally use, for example, the rate equations for the particle 
number density instead of the collision entropy density. With these other choices, the values of the different quantities 
are shifted somewhat due to the way that the initial conditions are extracted but there is no qualitative different in 
the result. Our present choice has the distinct advantage that we can explicitly compare the different processes using 



the collision entropy density rates. This will become clear when we show the results in Sect. VI. 



V. INITIAL CONDITIONS 

To start the evolution, we use the same initial conditions for the gluon plasma as before Jllj] based on HIJING result 
for Au+Au collision. The initial conditions for the quarks (antiquarks) are obtained by taking a ratio of 0.14 for the 
number of initial quark (antiquark) to the initial total number of partons as done in [jl3[-|T5|] . The initial conditions 
are shown in Table 1. One sees that the initial quark collision times are long compared to those of the gluons both at 
RHIC and LHC. Especially at RHIC, the quark collision time is exceedingly long and so these particles are essentially 
free streaming initially. Taking these numbers as guides to how fast each particle species is going to equilibrate, we 
can be sure already of a two-stage equilibration scenario [0 . 



Initial Conditions 





RHIC 


LHC 


t (fm/c) 


0.70 


0.50 


T (GeV) 


0.50 


0.74 


e 0g (GeV/fm 3 ) 


3.20 


40.00 


e , (GeV/fm 3 ) 


0.63 


7.83 


n 09 (fm~ 3 ) 


2.15 


18.00 


n q (fm~ 3 ) 


0.42 


3.53 


lo g 


0.08 


0.21 


h q 
60 g (fm/c) 


0.017 


0.044 


2.18 


0.73 


B q (fm/c) 


239.72 


30.92 



TABLE 1. Initial conditions for the evolution of a QCD plasma created in Au+Au collision at RHIC and at LHC 
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Using the standard initial picture of heavy ion collisions as before, our evolution is started when the momentum 
distribution in the central region of the collision becomes, for the first time, isotropic due to longitudinal cooling. The 
subsequent development is determined by the interactions Eq. (|ll|), ( p0| ) and (^l]). In the case of a pure gluon plasma 
|llf , it is clear that interactions bring the system towards equilibrium and not towards some free streaming final state 
which is a possible alternative as can be inferred from the analysis in Sect. III. That is the interactions dominate 
over the expansion. In the present situation, we will see that the same can certainly be said for the gluons and for 
the quarks at LHC but at RHIC, it is less clear for the latters. The equilibration time for quarks is at least several 
times longer than that of the gluons. 

Details for the procedure of the computation can be found in The values for the numerical parameters are 

the same and in addition, we use nf = 2.5 to take into account of the reduced phase space of strange quark. All 
time integrations are discretizcd and the rates are obtained at each time step necessary for forming the two pairs of 
equations Eq. ( p2f ) , ( |25| ) and ( prj| ) . One then solves the two equilibrium temperatures T eq g and T eq q from two 4th 
degree polynomials, one for each of the temperatures. From these solutions, 9 g and 9 q are obtained and everything is 
then fed back into the equations for the next time step. 



VI. EQUILIBRATION OF THE QCD PLASMA 



We show the results of our computation in this section. They show clearly the collision times 8 g and 9 q hold the 



keys to equilibration as have been analysed in Sect. III. We will see shortly that as a result of the disparity between 
their magnitudes at finite values of r, the equilibration of quarks and antiquarks lags behind that of the gluons both 
chemically and kinetically. We will also identify the dominant processes responsible for equilibration. They are not 
the commonly assumed elastic scattering processes as already mentioned in the introduction. 

When dealing with two particle species, one has several choices as to when should the evolution be stopped. We 
choose to do this when both the quark and the gluon temperature estimates drop to 200 MeV. For gluons, this 
estimate is obtained by the near equilibrium energy and number density expression 

e g = a 2 l g T g and n g = oi l g T g 3 , (35) 

which are valid when the fugacity l g is near 1.0 i.e. when the distribution functions can be approximated by 
f g (p,l g ,r) = l g fg(p,l g = 1 , r) . For quarks and antiquarks, we cannot do the same as l q has not time to rise 
above 0.5 so instead, the temperature is estimated from the same quantities in kinetic equilibrium but at small values 
of l q 

e q = 'iv q l q T q /n 2 and n q = v q l q T q /n 2 . (36) 




T (fm/c) T (fm/c) 

FIG. 1. The time-dependence of the estimated temperatures for quarks and for gluons and their fugacities at (a) LHC and (b) 
RHIC. The solid lines are the estimated temperatures T g (thick line) and T q . The dashed lines are the fugacities l g (thick line) 
and l q . Gluon chemical equilibration is much faster than that of the quarks. The curves are stopped when all the temperature 
estimates drop to 200 MeV. The vertical line indicates when the gluon temperature reaches this value. 
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These estimates are plotted in Fig. [jj The vertical line marks the point when the gluon temperature estimate (thick 
solid line) drops to 200 MeV. At this point, r ~ 6.25 frn/c, the fugacity (thick dashed line) is l g ~ 0.935 at LHC and 
is l g ~ 0.487 at r ~ 2.85 fm/c at RHIC. On the same plots, the quark temperature (solid line) drops at a slower rate 
and the fermionic fugacity (dashed line) is also increasing much slower given the less favourable initial conditions and 
initially much slower quark-antiquark pair creation than gluon multiplication rate. In the end, the fermions are not 
too well chemically equilibrated and in fact, are still quite far away from 1.0. This is especially bad at RHIC. We 
note that comparing to |15j-|lq], in our case, gluons chemically equilibrate faster but quarks are slower. 

Unlike chemical equilibration, kinetic equilibration has no simple indicators like the fugacities that can allow itself 
to be simply quantified. One has to, instead, use the anisotropy of momentum distribution as well as various reaction 
rates to get an idea of the degree of kinetic equilibration. The former can be deduced from the ratios of the longitudinal 
pressure and a third of the energy density to the transverse pressure, Pl/pt and e/3pr respectively. Whereas from 
the elastic scattering rates, one can deduce roughly how close the distribution functions are to their equilibrium forms 
by virtue of the fact that in local kinetic equilibrium, these rates are zero. The pressure ratios Pl/pt (solid line) and 
e/3pr (dashed line) are plotted in Fig. || (a) and (a') for gluons, (b) and (b ! ) for quarks and (c) and (c') for the total 
sum. These ratios are indeed approaching 1.0, the expected value after thermalization, but at different rates. Gluons 
are clearly equilibrating much faster than quarks which proceed rather slowly. 
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FIG. 2. The ratios of the longitudinal pressure (solid line) and a third of the energy density (dashed line) to the transverse 
pressure, Pl/pt and e/3pr respectively for (a) gluons, (b) quarks and (c) the total sum at LHC. Graphs (a'), (b') and (c') are 
the same at RHIC. 



To show that these behaviours, although slow, are indeed the signs of equilibration and that the plasma is not 
approaching some free streaming final states, we can work out what their behaviours should be in the latter case by 
taking the extreme and let 9i — > oo. From Eq. ([n]), as r — > oo, 
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6 -> 7TT e /4r J L ' ^ J ' 

where eo is the initial energy density and the above ratios are valid for both quarks and gluons in this extreme. 
Therefore in the free streaming case, the first ratio should approach zero and the second should approach 2/3. These 
are clearly not what we see in our plots. 




FIG. 3. The scaled products of the collective variables (a) energy density, (b) number density and (c) entropy density and 
their expected inverse time-dependence in equilibrium r 4/ ^ 3 , r and r respectively at LHC. Graphs (a'), (b') and (c') are the 
same at RHIC. The solid and dashed lines are for gluons and quarks respectively. The thick solid line in (c) and (c') is the 
scaled product of the total entropy density and r. 



To best get an idea of how close the distribution functions are to the equilibrium forms, the gg and qq or qq elastic 
scattering processes are ideal for this. These are shown in Fig. || and Fig. (b) for gluon and Fig. || and Fig. || (c) 
for quark. Note that the peaks of these collision entropy rates coincide with the corresponding mininum points of the 
pressure ratios. As expected, the rates maximize at maximum anisotropy in momentum distribution. They all rise 
rapidly from zero at To when the interactions are turned on. The subsequent return to zero or the approach of the 
distribution functions to their equilibrium forms are, however, much less rapid. They only do so progressively as can 
be deduced already from the pressure ratio plots. 

Having shown chemical and kinetic equilibrations separately, we present now the actual approach of the collective 
variables towards the equilibrium values. Since we are more interested in the behaviour of their time-dependence than 
their absolute magnitudes, we multiplied them by their expected time-dependence and scaled these by taking a guess 
at the corresponding asymptotic values from the tendency of the curves. The results are plotted in Fig. |[ They are 
£iT 4 ^ 3 /c S tTg{ 3 , riiT /n S iT S i and stT / s S iT S i in the figures (a) and (a'), (b) and (b') and (c) and (c') respectively. All 
these should be nearly constant with respect to time at large r. The solid lines are for gluons and the dashed ones are 
for quarks. They showed that the curves do behave in such a way for the eventual constant behaviour. This feature is 
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much clearer at LHC than at RHIC which only reconfirms the previously deduced result of faster equilibration at LHC 
than at RHIC. Note that for gluons, the quantities are approaching the corresponding asymptotic values from above, 
whereas for quarks, this approach is from below. This is because of the simple reason that there is a net conversion 
of gluons into quark-antiquark pairs via gg < — ► qq. The corresponding collision entropy density rate is negative as 
shown in Fig. || and Fig. ^ (c). We will see that this same interaction becomes dominant in the later part of the 
evolution later on when we compare the importance of the different processes. So gluons are losing energy, number 
and entropy to the fermions. This has to be so before the system as a whole can settle into complete equilibrium. 
The thick solid lines in Fig. ^ (c) and (c') show the scaled total entropy per unit area in the central region which give 
an idea of the state of the system as a whole. They show that although the entropy of the individual subsystem can 
decrease, the total value must increase in accordance with the second law of thermodynamics. 
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FIG. 4. The time-dependence of the collision time (a) for gluons g and (b) for quarks 9 q at LHC. Their values are compared 
in (c). t overtakes first 8 g and later 6 q also. Graphs (a'), (b') and (c') are the same at RHIC. In this case, r only has time to 
overtake 9 g but not 9 q . 



The figures discussed above show that the plasma is indeed approaching equilibrium and that interactions are fast 
enough to dominate over the Bjorken type one-dimensional scaling expansion. 

As we analysed in Sect. [II, thermalization is governed by the d^s. How fast this will proceed depends on their 
magnitudes and what is the actual final state depends on their time-dependent behaviours. For thermalization, the 
#i's must behave in such a way such that Xi — ► oo as r — > oo. That means they must grow less fast than r. In Fig. ^, 
we show these O^s as a function of r. Initially, 0j > r for both quarks and gluons, and 8 q starts off very large (see 
Table 1) but drops extremely rapidly back down to within hadronic timescales. The subsequent expected increase in 
time f§-||] is sufficiently slow for r to get past 6 g and 6 q at LHC, Fig. | (a) and (b) but at RHIC, Fig. | (b'), 6 q 
is still too large for r to overtake it before the temperature reaches 200 MeV. Nevertheless, the r-dependence is slow 
enough that Xi should go to infinity as r — > oo. 

We have mentioned in Sect. |l[ for the system to equilibrate as one, the target equilibrium temperatures T eq g 
and T eq q and also 8 g and 9 q must approach each other at large r. We strongly suspect that the convergence of the 
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temperatures will proceed in an oscillating fashion where the two curves intersect each other several times before 
the final convergence at very large r. We can see this in Fig. || (a) and (b). At LHC, the initial condition is more 
favourable for equilibration and so T eq g intersects T eq q twice already. This is not so at RHIC. In fact, all indications 
point to the fact that a plasma created at LHC will equilibrate better than one created at RHIC. By letting the 
plasma to continue its evolution and ignoring the deconfinement phase transition, we have seen that the collective 
variables like the gluon and quark energy densities, gluon entropy density etc. do show tendency to pass from below 
to above or vice versa, the corresponding equilibrium target values i.e. tendency to overshoot the equilibrium values 
and hence oscillation. As to the convergence of 0,'s, it is not so clear in Fig. |I] (c) and (c') 7 especially at RHIC in 
Fig. [| (c'). q is much too large in comparison with 9 g for any clear sign of convergence within the time available. 
On the other hand, at LHC, although there is still a large gap between the magnitudes, there is a clear tendency 
that the rate of increase of q with r is slowing down in Fig. [I| (b) while 6 g still increases at approximately the same 
rate. It is simply too early for the system to equilibrate as one. Even near the end, the quarks and gluons can only 
be considered as two linked subsystems approaching equilibrium at very different rates. Hence we have a two-stage 
equilibration. 
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FIG. 5. The time development of the equilibrium target gluon (solid line) and quark (dashed line) temperatures T eq g and 
T eq q respectively at (a) LHC and (b) RHIC. They should converge in an oscillating fashion at large r in order for the system 
to equilibrate as one towards a single temperature. The convergence is less good at RHIC than at LHC. 



Having shown that interactions can indeed dominate over the one-dimensional expansion of the parton gas in the 
central region of relativistic heavy ion collisions and hence bring the plasma into equilibrium. We can now look at the 
individual processes and compare their relative importance. These are the processes Eq. ©, © and ©. We have 
labelled their contributions to the gluon and quark collision entropy rate ds g /dr and ds q /dr by ds g i/dr, i = 1, . . .,4 
and ds q i/dT, i = 1, . ..,3 in the order that they appear in Eq. (19), ( |2C| ) and (^T|). Processes that give the same rate 
due to quark-antiquark symmetry are considered as the same process. Hence gq < — > gq and gq < — > gq give identical 
contribution to gluon and quark collision entropy density rate as dsgi/dr and ds^/dr respectively. Also we have 
combined fermion elastic scattering processes as one rate ds^/dr for convenience. These are shown in Fig. ^, Fig. [t], 
Fig. | and Fig. g. The elastic processes have a characteristic shape, i.e. an initial rapid rise to a peak at maximum 
anisotropy before returning to zero progressively. The sharper the peak, the quicker the kinetic equilibration (compare 
Fig. | and Fig. g (b), (d) and Fig. | and Fig. g (b), (c) and Fig. |). Note the negative rate of Fig. | and Fig. g (c) 
which is because there are net quark-antiquark pair creations from gluon-gluon annihilations and entropy decreases 
with the number of gluons as already mentioned in the previous paragraphs. We compare the different processes 
by plotting the ratio of the magnitude of each contribution to that of gluon multiplication for gluons in Fig. g and 
Fig. [7] (e) and the ratio of each rate to that of quark-antiquark creation for quarks (antiquarks) in Fig. || and Fig. g 
(d). In the (e) figures, gluon multiplication clearly dominates initially at r < 2 fm/c at LHC and r < 4 fm/c at 
RHIC since all three ratios in each plot are less than 1. After these times, qq creation becomes dominant (thick solid 
line) and rises to several times larger than gluon multiplication. The gg elastic scattering, on the other hand, tends 
to maintain a small, nearly constant ratio with gluon multiplication (solid line), which supports the claim made in 
That is, in a pure gluon plasma, gluon multiplication dominates over gg elastic scattering in driving the plasma 
towards equilibrium. This remains the case even when l g ~ 0.93 which shows that this dominance is not sensitive 
to the value of l g . The remaining ratio of quark-gluon scattering to gluon multiplication continues to rise but not as 



12 



rapidly as the first ratio. For quark entropy, Fig. || and Fig. ^ (d), both ratios of quark-gluon scattering (solid line) 
and fermion-fermion scatterings (dashed line) to gg < — > qq rate remain small during the time available although they 
are both on the rise. So for gluons, gluon multiplication dominates initially but is later overtaken by gg < — ► qq which 
continues to dominate over other elastic processes. For quarks (antiquarks), this same process dominates during the 
lifetime of the plasma. 

These behaviours can be understood in the following way. Gluon branching dominates initially over any other 
processes so long as gluons are not near equilibrium. Once they approach saturation (the l g estimates slow down 
their approach towards 1.0 in Fig. |l| (a) and (b) at about the times mentioned above), gluon-gluon annihilation to 
quark-antiquark takes over as the dominant one because the fermions are still far from full equilibration. Because of 
the latter reason, the other ratios involving quark or antiquark to gluon branching continue to rise. 




T (fm/c) 1 (fm/c) T (fm/c) 

FIG. 6. The time development of the different contributions to the total gluon collision entropy density rate at LHC. They 

are (a) gg < > ggg, (b) gg < > gg, (c) gg < > qq and (d) gq < > gq or gq < > gq. The curves of the elastic scattering 

processes in (b) and (d) have typical peaks at maximum anisotropy in momentum distributions. The ratios of the contribution 
(b) (thick line), (c) (solid line) and (d) (dashed line) to that of (a) are plotted in (e). This shows that first gluon multiplication 
dominates initially but is later overtaken by gluon annihilations into quark-antiquark pairs. 



So contrary to common assumption, inelastic processes are dominant in equilibration. This should have consequences 
in the perturbative calculations of transport coefficients or relaxation times pS6|~|3"8[] of system that are not subjected 
to external forces. These calculations are based essentially, up to the present, on elastic binary interactions. As we 
have seen, they are not the dominant processes in equilibration. 

To the surprising result of gluon multiplication dominates over elastic gluon-gluon scattering, we provide the 
following explanation. If one only looks at the scattering cross-sections, it is indeed true that gluon-gluon scattering 
has a larger value and gluon multiplication processes are down by a s for each extra gluon produced. The (n — 2) 
extra gluon production cross-section can be expressed in terms of the elastic scattering cross-section as 
double logarithmic approximation, 

o" 9ff ^(n-2)g oc <T gg -,gg [a s In 2 (s/s cu t)]™~ 4 (38) 

where s cu t is the cutoff for the mininum binary invariant (pi +Pj) 2 > s cu t of the 4-momenta of each gluon pair. In the 
present problem, s cu t = rn 2 D , the double logarithm is not large and certainly does not compensate for the smallness 
of ot s . However, as we have mentioned at the beginning, the collision term on the r.h.s. of Eq. (||) consists of the sum 
of the differences of the reactions in a QCD medium going forward and backward, so a large cross-section does not 
automatically imply dominance of the corresponding process in the approach to equilibrium. 
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FIG. 8. The time development of the different contributions to the total quark collision entropy density rate at LHC. They 

are (a) gg < > qq (b) gq < > gq or gq < > gq and (c) the sum of the contributions of all fermion elastic scattering processes 

qq < > qq, qq < > qq and qq < > qq. The ratios of the contribution (b) (solid line), (c) (dashed line) to that of (a) is plotted 

in (d). This shows that throughout the lifetime of the QCD plasma, gluon annihilations into quark-antiquark pairs dominates 
in the equilibration of the fermions. 

Similarly, gg < — > qq is not that different from gq < — > gq or gq < — > gq because the two matrix elements are related 
simply by a swapping of the Mandelstam variables. So why should the first dominates over the second? Except the 
different ways that the infrared divergences are cut off in the processes, the main reason is gg — > qq dominates over 
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the backward reaction qq — > gg due to the simple fact that there are less fermions than gluons present in the plasma. 
An extreme example of this phenomenon would be the forward and backward reaction balance out each other for all 
the elastic interactions as in a kinetically equilibrated plasma when only inelastic processes remain in the collision 
terms. In this extreme, all the ratios of elastic to inelastic collision entropy rate vanish. 

We can now return to the question of whether other inelastic processes such as gg < — ► qqg, gq < — ► gqg, gq < — ► gqg, 
gq < — ► qqq, gq < — ► qqq, qq < — ► qqg etc. should be included. Although they are non-leading compared to gg < — ► ggg 
and gg < — ► qq due to colour, they should be significant when one sizes them with the elastic processes in view of 
the cancellation between the forward and backward reactions. In |Tl[| , the question of the dominance of inelastic 
over elastic processes was raised. Here it is sufficient to include the two leading inelastic processes to show this 
explicitly. Had one included these other processes, then equilibration should be faster and one could end up with a 
more reasonable quark-antiquark content in the plasma. However, we are doubtful that the equilibration time can be 
reduced dramatically from what we have shown here. 
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As we argued in JiT[ , it is hard to perturb a parton system from thermal equilibrium without doing so chemically. 
Therefore inelastic processes are always active in the approach to equilibrium whereas the same is not true for elastic 
processes. From our figures, it can be seen that inelastic processes are not there only for chemical equilibration or 
for minor contributions to thermalization as is commonly assumed due to their possible higher powers in a s , they 
contribute even more significantly to equilibration than elastic processes. Changing the initial conditions will only 
vary the dominancy but not remove the dominance. 

Before closing, we would like to point out some differences of our results with that of PCM. In PCM, there appears 
to be no early momentaneous isotropic particle momentum distribution in either S+S or Au+Au collisions. The first 
time that there is approximate isotropy, it is already thermalization according to g. It was claimed that there was 
no further significant change in the total momentum distribution after r = 2.4 fm/c for Au+Au collision at RHIC. 
We assume that they mean the shape of the distribution with the exception of the slope which should continue to 
change due to cooling. However, when the total distribution is broken down into that of the parton components, the 
approximate isotropy or thermalization becomes less obvious. We have shown that thermalization in the strict sense 
is slow and isotropy of gluon momentum distribution can be argued to be approximate but that of the fermions is 
not so good. 

As to chemical equilibration, PCM shows little chance of that for the fermions. The corresponding fugacity estimates 
are approaching the "wrong direction" with increasing time. This is due to a net outflow of particles from the defined 
central region. The net flux of outgoing particles is arguably more important for fermions than for gluons because the 
formers have a larger mean free path. The result is the gluon (fermion) fraction of the particle composition rises (drops) 
with increasing time. Therefore even if there is no phase transition and the parton plasma is allowed to continue its 
one-dimensional expansion indefinitely, chemical equilibration will never be achieved. Then according to PCM, the 
expansion is slow enough for kinetic equilibration for all particle species but too fast for chemical equilibration of the 
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quarks and antiquarks. The boundary effect is too important and is affecting equilibration. In our case, this effect is 
not incorporated. Although equilibration is slow, full equilibration will be reached given sufficient time. 

We find it surprising that although the gluon fugacity estimate in PCM B overshoots and stays above or at 1.0 
nearly all the time except at the beginning, R g is still positive or an order of magnitude larger than R q + Rq when 
the fugacities of the latter are well below 1.0 and decreasing. One would expect rather gluon absorption or conversion 
into quark-antiquark should take a significant toll on the gluon production so that there should be a diminution of 
gluons. At least, this should be the case when local kinetic equilibrium has been or nearly been reached which PCM 
claimed to be so at the end of the program run but this is not the case in the plot of the production rate of the 
different particle species! This is counter-intuitive and opposite to what we have shown. 

To conclude, we have shown that inelastic processes dominate in the approach towards equilibrium. In particular, 
gluon branching is most important. Gluon-gluon annihilation into quark-antiquark becomes more important only 
when the gluons are near saturation and equilibrium. The lower power in a s of the gluon-gluon elastic scattering as 
compared to the inelastic gluon emission process is more than compensated for by the cancellation of the reaction going 
forward and backward. The recovery of isotropy in momentum distribution is slow and so is chemical equilibration. 
The latter is partly due to the small initial fugacities that we used. As an intrinsic feature of perturbative QCD, the 
quarks and antiquarks are lagging behind the gluons in equilibration and hence a two-stage equilibration scenario. 
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